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ABSTRACT 

Within a density-dependent relativistic mean-field model using in-medium meson-hadron cou- 
pling constants and meson masses, we explore effects of in-medium hyperon interactions on prop- 
erties of neutron stars. It is found that the hyperonic constituents in large-mass neutron stars 
can not be simply ruled out, while the recently measured mass of the millisecond pulsar J1614- 
2230 can constrain significantly the in-medium hyperon interactions. Moreover, effects of nuclear 
symmetry energy on hyperonization in neutron stars are also discussed. 

Subject headings: dense matter - equation of state - stars: neutron 



1. Introduction 

The in-medium hyperon interactions play an 
important role in determining properties of hy- 
pernuclei and the hyperonization in neutron stars. 
Conversely, observed properties of hypernuclei 
and neutron stars can be used to constrain the 
in-medium hyperon interactions. For instance, 
it has been shown that the properties of hyper- 
nuclei are indeed very useful in extracting in- 
medium hyperon potentials at subsaturation den- 
sities (Friedman & Gal 2007). In bulk matter, 
such as neutron stars, hyperons can be produced 
by virtue of strong interactions. They can actually 
become important constituents of neutron stars 
and thus have important effects on astrophysi- 



cal observations (for a review, see (Glendenning 
2001)). In fact, it is well known that the hy- 
peronization can reduce the maximum mass of 
neutron stars as much as 3/4M (Glendenning 
1985; Glendenning & Moszkowski 1991; Jiang 
2006). Currently, a number of phenomcnologi- 
cal models, considering only the minimum com- 
positions of nuclcons and leptons using inter- 
action parameters that are well calibrated by 
using terrestrial nuclear laboratory data, can 
produce the maximum mass of neutron stars 
around or below 2M Q (Danielewicz & Lacey 
2002; Piekarewicz 2007). Interestingly, several 
neutron stars with large masses around 2Mq have 
been observed (Nice et al. 2005, 2008; Ozel 2006) 
recently. In particular, the 2M@ pulsar J1614- 
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2230 was measured rather accurately through the 
Shapiro delay (Demorest et al. 2010). Since prop- 
erties of neutron stars are determined by the nu- 
clear equation of state (EOS) and the hyperoniza- 
tion reduces significantly the maximum mass of 
neutron stars, it has been stated that these ob- 
servations seem to rule out almost all currently 
proposed hyperon EOS. Recent evidence for this 
can also be found in the work of the Brueckner 
approach (Schulze & Rijken 2011). 

Though most of the hyperon EOSs give lighter 
neutron stars, there were actually a few endeavors 
in the past to stiffen the EOS either by invoking 
strong repulsions for hyperons or pushing upwards 
the onset density of hyperons, leading to heavy 
neutron stars involving hyperons (Hofmann et al. 
2001; Takatsuka et al. 2002; Stone et al. 2007; 
Dcxhcimer & Schramm 2008) . Recently, by virtue 
of nonlinear self-interacting terms involving a vec- 
tor meson with hidden strangeness, Bednarek et 
al. obtained a stiff hyperon EOS with which the 
large mass of the PSR J1614-2230 can be pro- 
duced (Bednarek et al. 2011). Usually, the SU(6) 
relations are imposed to constrain the meson- 
hyperon coupling constants. Such SU(6) rela- 
tions were recently reexamined by Weisscnborn 
et al. and an arbitrary breaking of such re- 
lations can also result in the stiffening of hy- 
peron EOS and the uplift of the maximum mass 
of neutron stars with the inclusion of hyper- 
ons (Weissenborn et al. 2011a). On the other 
hand, the quark deconfinement may occur in the 
medium as the spatial overlap of nucleons be- 
comes sufficient to dissolve the boundary of color 
singlets with the increase of density. Of course, 
the quantitative understanding of such a color 
deconfinement in cold medium is still model de- 
pendent. Typical effective QCD models include 
the NJL-like models, e.g., see (Klahn et al. 2007; 
Ippolito et al. 2008; Pagliara & Schaffner 2008; 
Bonanno & Sedrakian 2011) and widely used 
MIT bag models, e.g., see (Prakash et al. 1997; 
Alford et al. 2005; Weissenborn et al. 2011b), as 
well as Schwinger-Dyson approaches (Li and Luo et al. 
2011). With the Maxwell or Gibbs constructions 
for the hadron-quark phase coexistence, the result- 
ing quark EOS may give rise to two possible types 
of stars: strange stars which are totally made 
of absolutely stable strange matter (Glcndenning 
2000) and hybrid stars with a quark core and 



hadron out-layer. In order to be consistent with 
the recent observation of the 2Mq pulsar, the 
strong coupling and/or color superconductivity 
were shown to be necessary (Xu 2003; Alford et al. 
2005; Klahn et al. 2007; Ippolito et al. 2008; 
Pagliara & Schaffner 2008; Weissenborn et al. 
2011b; Bonanno & Sedrakian 2011). While the 
compositions of hybrid stars are rather model de- 
pendent, it is interesting to mention that Yasutake 
et al. suggested the hyperon suppression with the 
MIT model using a density-dependent bag con- 
stant (Yasutake et al. 2011). 

Noticing the recent investigations on the con- 
sequences of various quark EOS, we examine in 
this work the consistency of the hyperon EOS 
with the recent observation of the PSR J1614- 
2230. Despite the impressive progress made in 
recent decades in constraining the nuclear EOS 
using both astrophysical observations and nu- 
clear reaction data, see, e.g. (Youngblood et al. 
1999; Danielewicz & Lacey 2002; Li et al. 2008), 
many uncertainties still remain. The in-medium 
hyperon interactions are among the most un- 
certain ingredients of neutron star models. We 
shall thus seek in-medium hyperon interactions 
that can produce the observed maximum mass 
of neutron stars. In view of the fact that some 
microscopic theories, such as the Brueckner ap- 
proach, are still having difficulties to obtain the 
2Mq of hyperonized neutron stars, here we re- 
sort to the phenomenological models developed 
in refs. (Jiang et al. 2007a,b) to analyze effects of 
various in-medium hyperon interactions on prop- 
erties of neutron stars. 

2. Density-dependent relativistic mean- 
field models and parametrizations 

In order to study conveniently the in-medium 
interactions for hyperons, we seek for density- 
dependent relativistic models without nonlinear 
interactions. In our previous works (Jiang et al. 
2007a, b), we constructed density dependent rela- 
tivistic mean-field (RMF) models using in-medium 
hadron properties according to the Brown-Rho 
scaling due to the chiral symmetry restora- 
tion at high densities (Brown & Rho 1991; Song 
2001; Brown & Rho 2005; Brown et al. 2007). 
In these models, the symmetric part of the re- 
sulting equations of state around normal den- 
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sity is consistent with the data of nuclear giant 
monopole resonances (Youngblood et al. 1999) 
and at supra-normal densities it is constrained 
by the collective flow data from high energy 
heavy- ion reactions (Daniclcwicz & Lacey 2002), 
while the resulting density dependence of the 
symmetry energy at sub-saturation densities 
agrees with that extracted from the isospin dif- 
fusion data from intermediate energy heavy- 
ion reactions (Tsang et al. 2004; Chen et al. 2005; 
Li & Chen 2005). Our models with the chiral lim- 
its are soft at intermediate densities but stiff at 
high densities naturally, producing a heavy maxi- 
mum neutron star mass around 2Mq. It is inter- 
esting to see that the EOS extracted from the ce- 
lestial observations most recently features similar 
characters (Steiner et al. 2010). Apart from the 
usual studies on the minimum constituents in neu- 
tron stars with electrons, protons and neutrons, 
in this work we include the hyperonic degrees of 
freedom to study the in-medium interactions for 
hyperons with the constraints of the recent celes- 
tial observations. The model Lagrangian with the 
density-dependent parameters is written as 

C = i> B [i^-M* B +gl B a-gl Bllx ^ 

-g^Bl^K^B + -{d^ad^a - mfa 1 ) 

+im;V&o + £y + Ci, (1) 

where -0s, cr, w, and bo are the fields of the baryons, 
scalar, vector, and isovector- vector mesons, with 
their masses Mjlj , m* , to* , and to*, respectively. 
Fpu and B^ are the strength tensors of the uj 
and p mesons, respectively. The meson cou- 
pling constants and masses with asterisks denote 
the density dependence, given by the BR scal- 
ing (Jiang et al. 2007a,b). Li and Ly are the La- 
grangian for leptons and hyperons, respectively 
The parameters for strange mesons a* (i.e. /o, 
975MeV) and 4> (1020MeV) in C Y arc assumed to 
be density independent. 

The energy density and pressure in the RMF 
approximation read, respectively, 

c 1 *2 2 i 1 *2i2 i 1 2 i2 i *2 2 

£ = w + 2 TO p b o + 2 m ^o + if"-* ° 

w 



*2 2 i 1 *2i2 2 j.2 

P = w o + 2 TO p h o + 

1 *2 2 1 2 *2 

to fj m„,a — P 

2 <T rj (7 U ' 

where i runs over the species of baryons and lep- 
tons considered in neutron star matter, E* = 
^/k 2 + to* 2 with to* being the Fermion effective 
mass of species i, and is the rearrangement 
term, originating from the density-dependent pa- 
rameters, to preserve the thermodynamic consis- 
tency (Jiang et al. 2007a). 

The density dependence of parameters is de- 
scribed by the scaling functions that are the ra- 
tios of the in-medium parameters to those in 
the free space. For the nucleonic sector, we 
take the scaling functions for the coupling con- 
stants of scalar and vector mesons as ^aN(p) = 
(1 — yip/ po) I '(1 + Xip/po) with coefficients Xi and 
iji given in (Jiang et al. 2007a, b). For the hy- 
peronic sector, we will give the scaling form for 
coupling constants below. For hadron masses, the 
scaling function is given as $(p) = 1 — yp/po- 

For the neutron star matter with hyperoniza- 
tions, the chemical equilibrium is established on 
the weak interactions of baryons and leptons. We 
study chemically equilibrated and charge neutral 
matter including baryons (N, A, E, S) and leptons 
(e, p). Note that although the baryon chemical po- 
tential is modified by the rearrangement term, the 
chemical equilibrium is independent of this modi- 
fication. 

In the present work, the RMF parameter sets 
SLC and SLCd (Jiang et al. 2007b) that have no 
mass scalings for baryons are extended to in- 
clude the hyperonizations. The main difference 
between the SLC and SLCd is their prediction for 
the density dependence of nuclear symmetry en- 
ergy. More quantitatively, the slope parameter L 
at saturation density is L — 92.3 MeV and 61.5 
MeV for the SLC and SLCd, respectively. As- 
suming all analyses of terrestrial data are equally 
reliable, a conservative estimate puts L in the 
range of approximately 25 to 115 MeV and the 
symmetry energy at normal density E sym (po) to 
be between 26 and 34 MeV (Newton et al. 2011; 
Tsang et al. 2012). However, it is interesting to 
note that the majority of the analyses of ter- 
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restrial nuclear experiments scatter the L and 
Esym(po) values around 60 MeV and 30 MeV, re- 
spectively. Several recent analyses of astrophysical 
observations and phenomena including the mass- 
radius correlation (Steiner & Gandolfi 2012), the 
binding energy of neutron stars (Newton & Li 
2009), the frequencies of torsional crustal vi- 
brations (Steiner & Watts 2009; Gearheart et al. 

2011) and the r-model instability window (Wen et al. 

2012) , all consistently favor L values less than 
about 70 MeV. For instance, the latest analysis 
of neutron star observation puts L in a somewhat 
lower range of 36 to 55 MeV at 95% confidence 
level (Steiner & Gandolfi 2012). It is worth not- 
ing that the L value from the SLCd model is more 
consistent with the most stringent constraint on 
L so far, i.e. approximately 40 < L < 60 MeV, 
obtained by combining constraints from analyzing 
both terrestrial experiments and astrophysical ob- 
servations (Lattimer & Lim 2012). By comparing 
calculations with the SLC and SCLd parameter 
sets we shall examine effects of the symmetry en- 
ergy on the mass-radius correlation of neutron 
stars, though the SLCd model can describe the 
available observations more consistent with con- 
clusions based on other analyses in the litera- 
ture. The parametrization for hyperons is elab- 
orated in the following. The coupling of mesons 
with hyperons can generally be given in terms 
of the parameters X a y ^ X^y, and X p y, which 
are ratios of the meson coupling with hyperons 
to that with nucleons. Lack of strong constraints 
on these parameters (Weisscnborn et al. 2011a), 
a variety of choices of these parameters roughly 
varying from 0.2 to unity were used in prac- 
tical studies (Glendenning & Moszkowski 1991; 
Avancini & Menezes 2006). Considering that the 
nucleonic sector of our models respects chiral lim- 
its at high densities, we assume for hyperons two 
cases: the usual case (UC) that the hyperons 
have a similar medium effect to nucleons and the 
separable case (SC) that the meson-hyperon cou- 
pling constant is separated into density-dependent 
and density-independent parts regardless of the 
chiral limit constraint on the strange sector in 
hyperons. Nevertheless, the hypcron poten- 
tials (Millener et al. 1988; Hausmann & Weise 
1989; Fukuda et al. 1998) 

U[ N) = -30MeV = -UW, UL N) = -18MeV, 

(4) 



in nuclear matter at saturation density are used to 
preserve the relation between the vector and scalar 
meson coupling constants. Note that the repulsive 
£ hyperon potential is invoked here (Mares et al. 
1995; Noumi et al. 2002). For the strange mesons, 
we adopt the density-independent coupling con- 
stants in both cases for simplicity. The potentials 
for the A and S hyperons in S matter J7» = 
- = -40 MeV are used to obtain the 

coupling constants of the strange mesons (For the 
details, see (Jiang 2006; Schaffner & Gal 2000)). 
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Fig. 1. — (Color online) Coupling constants of the 
a and u mesons with the A hyperon as a function 
of density. 

The ratio parameters for coupling constants in 
the UC are assumed to be density-independent. In 
the SC, we consider following scaling functions for 
the meson-hyperon coupling constants that consist 
of two terms: 

1 2 
^a(p) = $tos(p) = ^^n(po) + 3^/v(p), 

, , 2 , , 1 
<Wp) = g$u,Jv(po) + 3*o,jv(p), (5) 

= (1 - faY)$<rN(Po) + Uy$vn{p), 

where p is the saturation density and f a y is an 
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adjustable constant. The scaling functions $ p s 
and $ p h for the p meson are taken the same 
as those of the ui meson. Note that $ w at(/3o) 
and $ CT iv(po) are just constants. The factors 
before $(po) and $(p) implied from constituent 
quark compositions play a role in averaging on 
the hadron level the coupling constant between 
the density-dependent part originating from the 
chiral limit and the density-independent part that 
is presumably attributed to the strange sector in 
hypcrons. The form in (5) produces the relation 
$iY = $i,/v at saturation density. Thus, we do not 
need to readjust the parameter gay(po) as f a y 
changes. For a few choices of / ct a, we plot the 
meson- A hyperon coupling constants in Fig. 1. It 
is seen that the larger the f a y , the smaller the g a Y 
at high densities. 

3. Results and discussions 

Since the parameters for the nuclconic sec- 
tor of the SLC and SLCd are clearly given 
in (Jiang et al. 2007b), we just list in table 1 the 
parameters for the hyperonic sector. We see in 
table 1 that in the same case all parameters but 
(?° s are the same for the models SLC and SLCd. 
This is because the unique difference between the 
models SLC and SLCd is that the latter has a 
softer symmetry energy than the former. Since 
the neutron star properties are rather insensitive 
to the coupling parameters of E and S hyperons 
owing to their small fractions in the core of neu- 
tron stars, for simplicity we take = fas = fa A 
for the SC in the calculation. For the similar rea- 
son, we prefer the choice X^a = A" W £ = X^s 
in the UC calculation, unless otherwise denoted. 
For the p meson in the UC, the usual relation 
X pT , = 2X pS = 2 is used. In both the SC and UC, 
the 3* A is zero. Note that all parameters used in 
this work for hypcrons meet the relation (4). 

As a well-known consequence, the emergence 
of the hyperon degree of freedom results in the 
softening of the EOS. While the models SLC and 
SLCd that were constructed based on the BR scal- 
ing, the softening turns out to be too appreciable 
at high densities to stabilize the neutron star in 
the UC with relatively small X^a- As shown in the 
upper panel of Fig. 2, this is related to the rapid 
decrease of the nucleon effective mass, correspond- 
ing to an increasingly large scalar field that pro- 



vides the attraction. As known before (Jiang et al. 
2007a), the vector coupling constant is decisive 
to generate a stiff EOS at high densities. Thus, 
the stiffening of the EOS can follow from increas- 
ing the parameter AT^a. With larger X u \, for in- 
stance, X^a = 0.9, the EOS is stiffened to recover 
the stability of neutron stars. For the SC, the EOS 
can be stiffened by increasing the parameter f a y, 
since the latter results in the decrease of the scalar 
coupling constant, as shown in Fig. 1. Generally, 
the EOS obtained with the SC is much stiffcr than 
that with the UC. Meanwhile, the accelerating de- 
crease of the baryon effective mass due to the in- 
clusion of hyperons can be greatly suppressed in 
the SC, as shown in the lower panel of Fig. 2. 
In Fig. 2, the A and S hyperon effective masses 
are also displayed. The upward shift of hyperon 
masses at high densities in the lower panel of Fig. 2 
is due to the decrease of the source term (namely, 
the hyperon density) of the strange mesons, also 
see below. The hyperon effective mass is much 
larger than the nucleon one. This would justify 
the use of the different in-medium interactions for 
hyperons and nucleons in the SC. 

The dropping of the baryon effective mass in 
chemically equilibrated matter associates tightly 
with the beginning density and fractions of hypcr- 
ons. In Fig. 3, we display the particle fractions as 
a function of density. We see that the fractions 
of hyperons are rather sensitive to the variation 
of the in-medium interactions. More noticeably, 
quite large differences between the left and right 
panels can be observed. For instance, in the UC 
(left panel), the interval between the A and S~ 
beginning densities depends sensitively on the ra- 
tios of coupling constants, although the relation 
(4) is always satisfied, while in the SC such a 
sensitivity does not exist. Moreover, the emer- 
gence of the E hyperon depends on the model 
interaction. The E hyperon does not appear in 
the SC. Remarkably, with the increase of density, 
the hyperon fractions in the SC tend downwards 
till to disappear after reaching the maximum, as 
shown in the right panel of Fig. 3. The occur- 
rence of this phenomenon is due to the fact that 
the vector meson-hyperon coupling constant <?* y 
has a weaker density dependence [see Eq.(5)] than 
the me son- nucleon coupling constant <7a;jv* 

At 

high densities, g* y exceeds g^ N , and so do the 
vector potentials. The chemical equilibrium thus 
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Table 1: Meson- hyperon coupling constants in various cases with models SLC and SLCd. In the UC, we 
tabulate the values with X u a = = = 0.9. The coupling constants in the medium are obtained 
as g* Y = g® Y ®iY(p)- In the UC, gf Y = 9iy(p = 0), while in the SC they are not equal. In the SC, the 
parameters listed here are free of the parameter j a y . For the couplings with the E hyperon, we take in the 
calculation the relations: = g° A and g° pTi = 2g° s . 
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Fig. 2. — (Color online) Baryon effective masses 
as a function of density in the UC and SC with 
the SLC. Three curves (solid, dashed and dotted) 
of the UC in the upper panel are calculated with 
Xu,z = 0.6 and X^a — 2/3,0.8 and 0.9, respec- 
tively. In the lower panel, three curves (solid, 
dashed and dotted) of the SC are obtained with 
faA = fas = 0.7, 0.8 and 0.9, respectively. 



makes the hyperon Fermi momenta and fractions 
lower. As an application in studying properties 
of neutron stars, this actually results in the ex- 
clusion of hyperons in the core of neutron stars 
and accordingly the re-stiffening of the EOS at 
high densities. In addition to our scheme, we 
note that there are other attempts to decrease 
the number of hyperons in neutron stars. For 
instance, different couplings of a new boson to hy- 
perons and nucleons were proposed to decrease 
the number of hyperons (Krivoruchenko et al. 
2009). Including the nonlinear self-interactions 
involving a vector meson with hidden strangeness, 
Bednarek ct al. found that the onset density 
of hyperons can be as high as 3/?o accompanied 
by smaller fractions of hyperons (Bednarek et al. 
2011). In (Takatsuka et al. 2002; Tsuruta et al. 
2009), however, hyperons were found to appear 
above 4po and thus played a rather limited role 
in the EOS of neutron star matter. In this work, 
we find that the onset density of hyperons in the 
SC can vary upwards within the region 2.2 — 3po 
as the ratio of vector meson coupling g^A/SuN 
increases from 0.2 to 0.8 (In table 1, this ratio is 
2/3), while the hyperon fractions get suppressed 
significantly at larger onset densities. A large 
rise of this density up to 4.5/50 can be obtained. 
However, the binding relation (4) should then be 
reduced to U A N) (p a ) = MeV. In this case, the 
hyperon fraction becomes so small that the hy- 
peronic constituents become unimportant. On 
the quark level, it is interesting to see that the 
mechanism of small numbers of strange quarks in 
hybrid stars was explored within a specific quark 
model (Buballa et al. 2004). 

We mention that the baryon fractions in neu- 
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Fig. 3. — (Color online) Particle fractions in the UC and SC with the model SLC as a function of density. 
In the UC (left panel), three curves for each particle, e.g., denoted by the number 1, 2, and 3, are calculated 
with X^a = 2/3, 0.8 and 0.9, respectively, while X u z = 0.9. The same number in the left panel denotes 
the results obtained with the same parameters. In the SC (right panel), three curves for each particle, for 
instance, in a rising order of the A hyperon appearance, are obtained with f a j\ = f a s = 0.7, 0.8 and 0.9, 
respectively. 



tron stars are also sensitive to the symmetry en- 
ergy. In (Jiang 2006), it was illustrated that the 
onset density for hyperons increases moderately 
with the softening of the symmetry energy. For 
the same reason, the onset density for hyperons 
with the SLCd is about 0.2p larger than those 
with the SLC. To save space herein, we do not 
display particle fractions with the SLCd in a fig- 
ure similar to Fig. 3. 

Although the emergence of hyperons is a cause 
for softening the nuclear EOS, the specific behav- 
ior relies indeed on the in-medium interactions, as 
clearly shown in Fig. 4. The softening persists at 
high densities for the UC. However, the softening 
is succeeded by a stiffening in the SC where hy- 
perons feel a different in-medium interaction from 
nuclcons. Looking back to the right panel of Fig. 3, 
we see that the stiffening of the EOS occurs with 
the suppression of hyperon fractions. As the hy- 
peron vanishes, the EOS returns to the normal 
EOS without hyperons. As shown in Fig. 4, the 



EOS evolves to become stiffer than the normal one 
with increasing density. The neutron star mat- 
ter thus transits to the normal isospin-asymmetric 
matter prior to the vanishing of hyperons. This 
eventually leaves a limited density window allow- 
ing the existence of hyperons in neutron stars. In- 
terestingly, we find that the influence of the hy- 
peronization in the SLCd model falls prominently, 
as compared to the SLC model. This is attributed 
to the larger onset density of hyperons with the 
SLCd due to the softening of the symmetry en- 
ergy, as mentioned above. 

In the UC, since the nucleon effective mass van- 
ishes at certain critical density, we need to con- 
sider the EOS beyond the critical density. Though 
the relationship between the chiral restoration and 
dcconfincmcnt occurrence is still discussed, it is 
nevertheless a convenient and usual way to ne- 
glect the distinction between them. Beyond the 
critical density, we thus adopt a quark matter 
EOS described by the MIT model with the ap- 
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Fig. 4. — EOS of isospin-asymmetric matter with 
models SLC and SLCd. Curves are obtained with 
parameters listed in Table 1. The arrow indicates 
the critical point to quark matter. The result in- 
cluding the muon but without hyperons is also de- 
picted for comparisons. 

propriate bag parameter [(179.5MeU) 4 ] to con- 
nect smoothly the hadron matter and quark mat- 
ter EOSs. The connection to the quark matter 
further softens the EOS, as shown in Fig. 4. 

We now turn to the consequences of hypcroniza- 
tions on properties of neutron stars. In partic- 
ular, it is interesting to see whether our model 
with hyperonization can give properties of neu- 
tron stars that are compatible with the recent 
observation of the millisecond pulsar J1614-2230. 
The mass of this pulsar was accurately determined 
to be 1.976 ± O.O4M (Demorest et al. 2010). It 
was concluded that such a large mass can rule 
out almost all currently proposed hyperon equa- 
tions of state. The mass-radius relation of neu- 
tron stars is obtained from solving the standard 
TOV equation with the nuclear EOS being spec- 
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Fig. 5. — (Color online) Mass-radius relation of 
neutron stars in the UC and SC with the SLC 
and SLCd. In the SC, two curves are obtained 
with f a Y = 0.8 and 0.9, respectively. The hatched 
areas give the probability distributions with la 
(red) and 2a (green) confidence limits for r p h 3> R 
summarized in (Steiner et al. 2010). 

ified above. For the low-density crust, we adopt 
the EOS in (Baym et al. 1971; Iida & Sato 1997). 
It is seen in Fig. 5 that all EOS's of the SC can 
produce neutron stars with the maximum mass 
around 2M Q and with the corresponding radius 
ranging from 9.3 to 9.5 km. Since the maximum 
mass of neutron stars is more sensitive to the 
EOS at high densities, it is understandable that 
the EOS including hyperonization at intermediate 
densities can lead to masses compatible with that 
of the pulsar J1614-2230. Thus, the EOS with 
hyperonization can not be simply ruled out. Of 
course, the hyperon fraction in neutron stars with 
the SC is rather limited, as compared to that with 
the UC. In addition, a subtle factor affecting the 
hyperonization is the density dependence of nu- 
clear symmetry energy. Due to the softening of 
the symmetry energy in the SLCd, the hyperoniza- 
tion in the SC becomes unimportant for the EOS 
of neutron star matter, see Fig. 4, which cquiva- 
lcntly expels most hyperons in neutron stars. Con- 
sequently, the mass-radius relation with the SLCd 
is almost independent of the parameter ] a y in the 
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SC, as shown in Fig. 5. For the UC, it looks un- 
doubtedly that the EOS is ruled out by the re- 
cent observation (Dcmorest et al. 2010), since the 
maximum mass of neutron stars in this case just 
sprints to 1.4M@ with a much compacter size than 
the canonical neutron star without hyperons. 

It is worth adding here some discussions about 
the influence of the symmetry energy on the 
radii of neutron stars. It is now well estab- 
lished that the maximum mass of neutron stars 
is dominated by the high-density behavior of the 
EOS, while the radius is primarily determined 
by the slope of the symmetry energy at inter- 
mediate densities (1 -3po) (Lattimer & Prakash 
2001, 2004; Stciner ct al. 2005; Xu et al. 2009). 
The large difference between the radii of low-mass 
neutron stars obtained with the SLC and SLCd 
can be attributed to the difference in the slope 
parameter L. As discussed in detail earlier in 
ref. (Fattoyev & Piekarewicz 2010), the central 
density of an 18 - 20 km star is near po, and the 
crust ends approximately at (1/3 — 1/2) pq. In 
this density range, the pressure is dominated by 
the symmetry energy and not the incompress- 
ibility of symmetric nuclear matter. Our re- 
sults shown in Fig. 5 are consistent with those 
in ref. (Fattoyev & Piekarewicz 2010). For the 
1AM@ neutron stars in the SC, we see however 
that the inclusion of hyperons reduces moderately 
the range of the neutron star radius. For instance, 
the radius ranges from 11.4km with the SLCd to 
12.3km with the SLC for f a y = 0.9, which are 
well situated in the domain extracted by Steiner 
et al. (Steiner et al. 2010). With the emergence 
of hyperons, the sensitivity of the star radius to 
the symmetry energy is reduced clearly. This is 
because the A-hyperon, being the dominant com- 
ponent of hyperons, is an isospin scalar. The sup- 
pression of the isovector potential g* p b$ due to the 
appearance of A hyperons (Jiang 2006) results in 
the reduction of the pressure of asymmetric mat- 
ter and thus the star radius, while the magnitude 
of the suppression depends on the specific values 
of the symmetry energy and hypcron fraction. As 
a result, the less sensitivity of the star radius to 
differences in the symmetry energy is observed in 
calculations with both the SC and UC. In Fig. 5, 
we also include the constraints of mass-radius tra- 
jectories for r p h 3> R obtained by Steiner et al. 
in (Steiner ct al. 2010). Our results in the SC are 



either within (for SLCd) or not very far off (for 
SLC) the constrained region, though the inclusion 
of hyperons seems to tilt the vertical trajectories. 
For low-mass neutron stars, the radii with the SLC 
are predicted to go beyond the optimal region ex- 
tracted very recently in (Steiner et al. 2012), un- 
less some other scenarios that allow a loose exten- 
sion of the radii are invoked to extract the radius 
constraints (Suleimanov et al. 2011; Zhang et al. 
2007). However, the maximum mass of neutron 
stars is almost independent of the slope L and 
the radii of low-mass neutron stars. In our cases, 
the maximum mass is only reduced by about 1% 
by softening the symmetry energy from the model 
SLC to SLCd. On the other hand, it is known that 
the measurements of the neutron star radii are far 
less precise than the mass measurements, see, e.g., 
refs. (Suleimanov et al. 2011; Zhang et al. 2007) 
and references therein. We note that an option of 
Tph = R can cause a visible slanting of the vertical 
trajectories (Steiner et al. 2010). The agreement 
of our results can thus be better with the con- 
straints obtained for r pc = R. 

We note that there had been a few endeavors 
in the past to involve the hyperons in heavy neu- 
tron stars either by invoking strong repulsions for 
hyperons or pushing upwards the onset density of 
hyperons(Hofmann ct al. 2001; Takatsuka et al. 
2002; Stone ct al. 2007; Dexheimer k Schramm 
2008; Bednarek et al. 2011; Wcisscnborn et al. 
2011a). The main purpose of our work is to con- 
strain the in-medium hyperon potentials using the 
2 solar mass constraint of neutron stars. This in- 
deed requires strong repulsions for hyperons. Due 
to different density-dependencies of nucleonic and 
hyperonic potentials in the SC, the hyperonic vec- 
tor potential exceeds the nucleonic one, leading to 
a very significant suppression of hyperons at high 
densities. As a result, the hyperons would exist in 
a shell of neutron star core even with a small L 
value. 

Besides the effect of hypcronization on the 
static properties, another consequence of the hy- 
pcronization is on the thermal evolution of neutron 
stars. For the SC in SLC, we see from Fig. 5 that A 
hyperons start to appear above 2.5po- At this cen- 
tral density the neutron star mass is about 1.2Mq. 
For most observed neutron stars that have larger 
masses, it appears that the direct Urea (DU) pro- 
cess with nucleons (Lattimer et al. 1991) and/or 



9 



hypcrons (Prakash ct al. 1992) would occur since 
the proton fraction in the neutron star matter 
with hyperonization can be in excess of the DU 
threshold (14%) with the SC. According to the 
thermal evolution of observed neutron stars an- 
alyzed in (Page et al. 2004; Yakovlev & Pethick 
2004), the fast cooling with the DU processes 
seemed to be excluded in most neutron stars ex- 
cept the massive ones. Slower cooling is possible 
when the neutrino emissivity can be suppressed 
by the superfluity of constituent particles like nu- 
cleons and hyperons when the temperature falls 
below the critical temperature. Page et al. set a 
stringent requirement for the critical temperature 
of neutron superfluidity without which enhanced 
cooling from DU processes may be needed in at 
least half of the observed young cooling neutron 
stars (Page ct al. 2009). While the DU cooling 
involving nucleons only was regarded to be too 
fast, it was pointed out by Tsuruta et al. that 
the DU cooling with hyperons in neutron stars 
can be compatible with the observations, provided 
that the hyperon superfluidity is appropriately ac- 
counted for (Tsuruta et al. 2009). In our models 
with the SC, the hyperonization can thus be favor- 
able to make up the potential incompatibility in 
the thermal evolution by considering the hyperon 
superfluidity. On the other hand, the threshold 
mass of neutron stars which allow the DU process 
increases with the softening of the symmetry en- 
ergy due to which the proton fraction exceeds the 
threshold value at larger densities. For instance, 
this mass is around 1.3M Q with the SC of the 
model SLCd, and the threshold density for the 
DU process is about 4po, which is coincidently the 
onset density of hyperons in (Tsuruta et al. 2009). 

Finally, we discuss some details concerning 
the in-medium interactions for hypcrons. Con- 
strained by the large mass of observed pulsars, 
it is favorable for us to select the in-medium 
interactions for hyperons in the SC. It is how- 
ever interesting to find that the single-particle 
potentials for hyperons in the UC and SC al- 
most overlap in a large density domain ranging 
from zero density to intermediate density. The 
significant departure appears only at high densi- 
ties (roughly > 4po)- Especially at lower densi- 
ties, the SC and UC descriptions give rather lim- 
ited differences in single-particle properties. This 
thus indicates that without compromising the suc- 



cess in describing properties of hypernuclei, one 
can constrain significantly the high-density hy- 
peron interactions with the large mass of neu- 
tron stars. In addition, we have noticed that 
many efforts using various quark models with the 
postulate of strong interactions and/or color su- 
perconductivity can also lead to large-mass hy- 
brid stars (Alford et al. 2005; Klahn et al. 2007; 
Ippolito et al. 2008; Weissenborn et al. 2011b; 
Bonanno & Sedrakian 2011) with stiff EOS's at 
high densities. Our models respecting chiral limits 
possess the similarly stiff EOS at high densities. 
On the other hand, neutron stars may be com- 
posed of more complicated constituents including 
quarks and meson condensates. It is useful to con- 
sider these non-baryonic degrees of freedom and 
their interplay with hyperons to constrain more 
quantitatively the in-medium interactions for hy- 
perons. However, this is beyond the scope of the 
present work. 

4. Summary 

The density-dependent relativistic mean-field 
model, which was constructed to respect the chi- 
ral symmetry restoration at high densities in terms 
of in-medium hadron properties according to the 
Brown-Rho scaling, is extended to include the hy- 
peronization in isospin-asymmctric matter. We 
examined effects of the in-medium hyperon inter- 
actions on properties of neutron stars. It is found 
that the maximum mass of neutron stars can con- 
strain significantly the in-medium hyperon inter- 
actions. In particular, assuming two categories of 
in-medium interactions for hyperons, we investi- 
gated their distinct roles in hyperonizations and 
properties of neutron stars. With different in- 
mcdium interactions for hyperons and nucleons 
(i.e., the SC case), a maximum neutron star mass 
of 2Mq can be obtained with the model where 
the nuclconic EOS is consistent with the terrestrial 
nuclear laboratory data. The result in the SC is 
compatible with recent observation on the mass of 
the millisecond pulsar J1614-2230. In this scheme, 
the number of hyperons is limited in neutron stars. 
Interestingly, we also found that the softening of 
the symmetry energy can play an important role 
in further reducing hyperons in neutron stars. On 
the contrary, the scheme that adopts similar in- 
medium interactions for both hyperons and nu- 
cleons (i.e., the UC case) can give a maximum 
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neutron star mass of only about 1.4M Q . Never- 
theless, the hypcronization in both cases reduces 
clearly the sensitivity of the neutron star radius to 
the difference in the nuclear symmetry energy. 
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